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We have investigated polycrystalline samples of the zigzag chain system BaTb 204 with a combination of 
magnetic susceptibility, heat capacity, neutron powder diffraction, and muon spin relaxation measurements. 
Despite the onset of Tb®"*" short-range antiferromagnetic correlations at |0cw| = 18.5 K and a very large 
effective moment, our combined measurements indicate that BaTb 2 04 remains paramagnetic down to 0.095 K. 
The magnetic properties of this material show striking similarities to the pyrochlore antiferromagnet Tb 2 Ti 207 , 
and therefore we propose that BaTb 204 is a new large moment spin liquid candidate. 

PACS numbers: 75.40.-s, 75.25.-j, 75.30.Cr 


Spin liquids are exotic ground states of frustrated magnets 
in which local moments are highly correlated but still fluc¬ 
tuate strongly down to zero temperature H]. In principle, the 
fluctuations of a spin liquid can be quantum or classical in na¬ 
ture. Several types of spin liquids have been proposed theoret¬ 
ically, including Anderson’s resonating valence bond stately, 
spin ice in |3 and others characterized by either gapped or 
gapless low-energy excitations Q. The experimental search 
for new spin liquid candidates is an ongoing area of inter¬ 
est since this state of matter remains largely unexplored in 
the laboratory. Some well-known examples of quantum spin 
liquid candidates include ZnCu3(OH)6Cl2(herbertsmithite)||6l 
111 , BaCu3V208(0H)2(vesignieite)|[Sl, and Ba3NiSb209ll9l, 
while their classical counterparts are the pyrochlore magnets 
Ho2Ti207 and Dy 2 Ti 2 O 7 llIl[l 0 llII]. 

Magnetic systems characterized by frustration and low- 
dimensionality have proven to be useful starting points in 
the quest for uncovering additional spin liquids, but partic¬ 
ular complications have severely limited the number of vi¬ 
able candidates. For example, although one-dimensional (ID) 
magnets and two-dimensional Heisenberg systems are not ex¬ 
pected to order due to the Mermin-Wagner theorem ifT^ . most 
real low-dimensional materials are governed by weak ex¬ 
change interactions in the other spatial dimensions and there¬ 
fore exhibit magnetic order at low temperatures. Furthermore, 
while magnetic frustration can drastically suppress leading 
terms in the Hamiltonian that would be responsible for mag¬ 
netic order, there are instances in which the sub-leading terms 
can still drive the system to an ordered ground state 1 ( 131 . To 
overcome these obstacles and And new spin liquid candidates, 
it is important to perform detailed studies on a wide variety of 
frustrated, low-dimensional magnets. 

The family of materials AR2O4 (A = Ba, Sr; R = rare 
earth) lfT4l - [T6l satisfy the two criteria described above. Two 
crystallographically-inequivalent R sites independently form 
two different types of zigzag chains running along the c-axis, 
and therefore quasi-ID magnetic behavior may be expected. 


Bulk characterization studies have also shown that most mem¬ 
bers of the family have dominant antiferromagnetic (ARM) 
exchange interactions and relatively large frustration indices. 
Geometric frustration can arise in this structure type if the J2 
exchange interactions are ARM and the Ji couplings are of 
comparable strength. 

The general trend for the AR2O4 family is that a small 
(large) ionic radius ifTTl r for the R atom shortens (lengthens) 
the Ji bonds and induces Neel (double Neel) order on the 
chains. Schematics of these two magnetic ground states are 
shown in Fig. [TJa). The Neel and double Neel long-range 
ordered states are consistent with predictions for the large 
ARM Ji and J2 limits of the classical Ising J1-J2 chain ifTSl . 
and have been realized in SrYb204(r = 0.87 A) l(T9l and 
BaNd204(r = 0.98 A) !^ respectively. By tuning the ionic 
radii of the R atoms with different rare earths, the intermedi¬ 
ate regime between these two limits has also been explored 
in detail. A variety of magnetic ground states have been 
observed, including coexisting long-range Neel and short- 
range double Neel order in SrEr204(r = 0.94 A)||2T1 1^ . 
coexisting short-range Neel and short-range double Neel or¬ 
der in SrHo204(r = 0.90 A) ll23l425l . incommensurate mag¬ 
netic order in SrTb204(r = 0.92 A) ll26l . short-range mag¬ 
netic order in SrDy204(r = 0.91 A) at 0.05 K|l22l ED, 
and no magnetic ordering of any kind down to 0.065 K in 
SrTm204(r = 0.88 A) ll29l . It is also interesting to note that 
any long-range order observed in this family has been found to 
arise from only one rare earth site (except for SrYb 904 l(T^ ). 
and therefore the two types of zigzag chains are often charac¬ 
terized by different magnetic ground states. Since the Ji and 
J2 bond lengths are essentially equal for the two chain types, 
this behavior may be a consequence of the inequivalent, dis¬ 
torted oxygen octahedral local environments of the rare earths 
comprising each chain. 

In this Letter, we investigate polycrystalline BaTb204 
(r = 0.92 A) with a combination of magnetic susceptibility, 
heat capacity, neutron diffraction, and muon spin relaxation 
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FIG. 1: (a) Schematics of the zigzag chain showing both the Neel and double Neel ground states that arise from the classical J 1 -J 2 Ising 

model in the strong AFM Ji and J 2 limits respectively, (b) Magnetic susceptibility of BaTb 204 , showing no evidence for long-range order 
down to 2 K. (c) Cp data measured between 0 and 12 T are shown on a double logarithmic scale with no indication of magnetic order. The low 
temperature maximum at T* = 1.5 K in zero field is likely due to a low-lying crystal field level. 


measurements. The magnetic species in this system is Tb^+, 
which is a non-Kramers ion with a large angular momentum 
J — 6. Despite the onset of antiferromagnetic correlations 
at the Curie-Weiss temperature 9cw = -18.5 K, there is no 
evidence for long-range magnetic ordering or spin freezing 
in any of the measurements down to 0.095 K. However, neu¬ 
tron diffraction reveals that incredibly short-ranged magnetic 
correlations exist between J 2 bonds. These findings provide 
strong evidence that BaTb 204 is a new large moment spin liq¬ 
uid candidate. 

Polycrystalline BaTb 204 samples were synthesized by a 
standard solid-state reaction method from high-purity starting 
materials of BaCOa and TbrOy. First, Tb 203 was obtained by 
reducing TbrOy in Ar (4% H 2 ). Next, a stoichiometric mix¬ 
ture of BaCOa and TbaOa (with 10% excess BaCOa) were 
ground, pressed into pellets, and then sintered in Ar (4% H 2 ) 
at 1150°C for 8 hours. The final product was confirmed to be 
single phase by laboratory x-ray powder diffraction. 

The magnetic susceptibility of polycrystalline BaTb 204 
was measured in an applied field fioH = 0.1 T under the 
zero-field-cooled condition using a Quantum Design Mag¬ 
netic Properties Measurement System. The data are presented 
in Fig.j^b), plotted as 1/%, i.e. ^qH/M, vs. T, and the results 
are in good agreement with previous work lfTSll . The high tem¬ 
perature data is well-described by a Curie-Weiss law, with a 
fit between 40 and 100 K yielding 6cw = -18.50(3) K and an 
effective moment /ie// = 9.95(1) fiB- The effective moment 
is close to the expected value of 9.72 for Tb^+. Despite 
the onset of AFM correlations around 20 K, there is no ev¬ 
idence for long-range magnetic order from the susceptibility 
measurements down to 2 K. 

The specific heat {Cp) below 2 K was measured with a 
home-built probe based on the adiabatic heat-pulse technique 
in a He-3/He-4 dilution refrigerator from Oxford Instruments. 
Data at higher temperatures were taken in a Quantum De¬ 
sign Physical Properties Measurement System equipped with 
a 12 T superconducting magnet. Fig.[TJc) shows Cp data in se¬ 
lected applied fields. No evidence for magnetic order is found. 


while te observed field dependence is likely indicative of crys¬ 
tal field splitting that changes with ^qH. 

The zero field Cp data show a sharp upturn for T < 0.3 K 
that can be attributed to nuclear Schottky contributions of Tb 
nuclei, as discussed for other Tb compounds ll^ 1311 . A broad 
maximum is also observed at T* = 1.5 K, which likely cor¬ 
responds to a low-lying crystal field level. Integrating CpjT 
from 0.5-6 K over the broad peak(assuming a negligible lattice 
contribution in this range|| 2 Tl) yields an entropy of just over 
i?ln(3/2)/mol-Tb. This finding could imply a doublet ground 
state for the Tb^"'' ions with a small energy gap to a low-lying 
singlet excited state. However, we note that Cp measurements 
on the AR 2 O 4 family provide limited quantitative information 
on the crystal field level schemes, due to the monoclinic site 
symmetries and the inequivalence of the local environments 
for the two different R sites. The build-up of short-range mag¬ 
netic correlations with decreasing T can also complicate anal¬ 
ysis of the magnetic entropy extracted from Cp data lf32l . 

Neutron powder diffraction was performed with 5 g of 
polycrystalline BaTb 204 between 0.3-100 K at Oak Ridge 
National Laboratory using the HB-2A powder diffractome¬ 
ter of the High Flux Isotope Reactor with a collimation of 
12 '-open- 6 '. Data with a neutron wavelength of 2.41 A is de¬ 
picted in Fig. I^a) and (d) with T = 0.3 K and 10 K respec¬ 
tively. Successful Rietveld refinements were performed using 
FullProf||33] with the known room temperature space group 
Pwaw lfTSll . indicating that there are no structural phase transi¬ 
tions down to 0.3 K. The lattice constants at 0.3 K refined as 
a = 10.423(1) A, 6 = 12.178(1) A, and c = 3.497(1) A. 

No evidence was found for long-range order in the diffrac¬ 
tion data of BaTb 204 . However, magnetic diffuse scattering 
was observed instead in both the 0.3 and 10 K datasets. This 
contribution is modeled as background in Fig. |^a) and (d), 
and most clearly seen in Fig.|^b) and (e). The Q-dependence 
of the diffuse scattering remains almost unchanged up to 10 K. 
These findings indicate that there are significant magnetic cor¬ 
relations that persist well above the onset of any possible long- 
range order. We note that the combined diffraction and sus- 
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FIG. 2: (a) HB-2A neutron diffraction data with A = 2.41 A at 0.3 K for polycrystalline BaTb 204 . The solid curve is a fit generated from a 
structural Rietveld refinement using the space group Pnam. (b) An enlarged version of the data shown in panel (a), emphasizing the magnetic 
diffuse scattering, (c) A 0.3-100 K difference plot, with the intensity normalized by the Tb^"*" magnetic form factor squared. The dashed and 
solid curves are fits to models incorporating only Tb®"*' magnetic correlations between Ji and J 2 bonds respectively, (d)-(f) Similar plots to 
those shown in panels (a)-(c), but with T = 10 K. 


ceptibility data rule out a well-isolated crystal field singlet 
ground state as found for Tm9Ti907 ll34l and possibly applica¬ 
ble to SrTm904 ll^ . since the signatures for such a scenario 
are a constant low-T magnetic susceptibility and the absence 
of any elastic magnetic scattering, in contrast to observations. 

0.3-100 K and 10-100 K difference plots are shown in 
Fig. I^c) and (f), with the data normalized by the Tb^+ form 
factor squared. Similar oscillatory scattering patterns are 
clearly evident in both datasets, therefore no drastic change is 
found in the magnetic correlations through T*. This observa¬ 
tion confirms that the Cp anomaly at T* is not associated with 
any form of magnetic ordering, but can likely be attributed to 
a low-lying crystal field level. The intensity of the difference 
plots is well-described by the function; 

( 1 ) 

ij 

This equation represents the expected polycrystalline re¬ 
sponse for a local magnetic structure with the spins at sites 
i and j correlated over distances of dij only. Antiferromag¬ 
netic (ferromagnetic) correlations are inferred from Aij < 0 
(Aij > 0). In BaTb 204 , the 0.3 K refinement discussed above 
yields a Ji bond length of 3.50 A, while the two types of 
zigzag chains are found to have inequivalent J 2 bond lengths 
of 3.59 A and 3.60 A. The best fits of the BaTb 204 data, 
shown by the solid curves in Fig. |^c) and (f), require only 
one term with A < 0 and d = 3.58(4) A at 0.3 K (3.61(3) A at 


10 K). Fits with dij fixed to 3.50 A are also shown for compar¬ 
ison purposes by the dashed curves, but they do simulate the 
data as accurately. Therefore, the analysis described above is 
consistent with AFM correlations extending only between J 2 
bonds, with d representing an average bond length between 
the two types of chains. Furthermore, the incredibly short- 
range nature of the correlations is reminiscent of the magnetic 
behavior in the large moment spin liquid candidate Tb 2 Ti 207 , 
where the Tb^+ spins were found to be correlated over a sin¬ 
gle tetrahedron onlv OSl . 

The neutron diffraction results discussed above do not al¬ 
low one to determine whether the spins are static or dynamic 
in the magnetic ground state. For example, diffuse elastic 
magnetic scattering was observed for all three pyrochlore sys¬ 
tems Y9Mo907 ll^ IJTI . Tb9Mo907 ll^ . and Tb9Ti907 ll^ . 
but the molybdates have been characterized as spin glasses 
while the titanate is best described as a spin liquid. Muon spin 
relaxation (/rSR) measurements proved to be instrumental in 
unambiguously determining the nature of the ground states in 
these cases. More specifically, the /iSR spectra for the two 
scenarios evolve quite differently in a longitudinal magnetic 
field confi guration BOl . 

Therefore, to better understand the magnetic ground state 
of BaTb 204 , pSR was performed at the EMU (10-300 K) and 
MuSR (0.095-4 K) spectrometers in longitudinal field (LF) 
geometry at the ISIS Pulsed Neutron and Muon Source. Zero 
field data were collected at 0.095, 0.5, 2, and 4 K, and no fea- 
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FIG. 3: (a) /iSR spectra collected in selected LFs at 0.095 K. Very 
little decoupling of the zero field relaxation is observed up to 0.2 T, 
which suggests that its origin is a dynamic mechanism, (b) T- 
dependence of the relaxation rate Ai, which saturates below 1 K. 
This T -dependence is consistent with dynamic behavior of the spins 
down to the lowest temperatures investigated. 

tures characteristic of long-range order or spin freezing were 
observed. Data were also collected in several different LFs at 
0.095 K. Some selected spectra are shown in Fig. [^a), and 
they fit well to the expression; 

A{t) = (2) 

where the two terms represent muons that stop in the sample 
and the sample holder respectively. The single sample com¬ 
ponent suggests that both types of chains host the same mag¬ 
netic ground state, in contrast to the picture that has emerged 
for most other AR 2 O 4 systems. This is consistent with the 
diffraction response that can be effectively modeled with only 
one local magnetic structure. 

One can estimate the local field at the muon site in 
the sample, assuming it is static, from the zero field /rSR 
data by using the approximation; Bioc = Ai/ 7 ^, where 
7 ^ = 135.5 X 27r MFlz/T is the muon gyromagnetic ratio. 
In the present case, this analysis yields Bioc = 0.0125 T for 
the lowest temperatures studied. Typically, if the observed 
muon spin relaxation in zero field is from a static mechanism, 
a LF up to one order of magnitude larger than the local field 


should be enough to completely decouple the relaxing signal. 
Fig. I^a) shows that the relaxation from the sample holder 
contribution is easily decoupled and therefore Xb =0 for the 
LF spectra presented in this work. However, inspection of 
Fig. [^a) shows that the opposite is true for the relaxation as¬ 
sociated with the sample contribution. In fact, the spectrum 
with hqHlf = 0.2 T shows only modest changes from the 
zero field case. This behavior is consistent with the muon 
spin relaxation arising from a dynamic origin and certainly 
not characteristic of a spin glass ground state. 

fiSR spectra were also collected with iiqHlf = 0.02 T at a 
variety of different temperatures to allow for an accurate mea¬ 
surement of Ai = 1/Ti. Fig.|^b) depicts the T-dependence 
of 1/Ti. No sharp peak is observed as would be expected 
in the case of spin freezing lWl . Instead, 1/Ti saturates be¬ 
low ~ 1 K, with a value roughly five times greater than found 
for Tb9Ti907 llT5ll . A single exponential muon spin relaxation 
from the sample is also observed at all temperatures. These 
findings are indicative of large, rapidly-fluctuating internal 
magnetic fields, as expected for a paramagnet BOl , and they 
are in good agreement with the decoupling scheme observed 
in the LF scan at 0.095 K. 

The absence of magnetic transitions in the bulk measure¬ 
ments, the presence of low-temperature short-range spin cor¬ 
relations from neutron diffraction, and the lack of any static 
magnetism detected by /rSR down to 0.095 K provide strong 
evidence for a large moment spin liquid ground state in 
BaTb 204 . This exotic state likely arises from incredibly bal¬ 
anced Ji and J 2 intrachain exchange interactions, which lead 
to a high degree of frustration. Single crystal inelastic neutron 
scattering and theoretical calculations are needed to help de¬ 
termine the crystal field level schemes and the evolution of the 
magnetic Hamiltonian throughout the AR2O4 series, as such 
work will lead to a deeper understanding of the origins of the 
exotic and varied magnetism in this family of materials. 
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